Simple, empirical order parameter for a first order quantum phase transition in 

atomic nuclei 
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A simple, empirical signature of a first order phase transition in atomic nuclei is presented, the 
ratio of the energy of the 6 + level of the ground state band to the energy of the first excited + 
state. This ratio provides an effective order parameter which is not only easy to measure, but also 
distinguishes between first and second order phase transitions and takes on a special value in the 
critical region. Data in the Nd-Dy region show these characteristics. In addition, a repeating degen- 
eracy between alternate yrast states and successive excited + states is found to correspond closely 
to the line of a first order phase transition in the framework of the Interacting Boson Approximation 
(IBA) model in the large N limit, pointing to a possible underlying symmetry in the critical region. 



The study of structural evolution in atomic nuclei has 
witnessed significant developments in recent years. One 
of the most important has been the discovery of empiri- 
cal evidence [T|, |2| for quantum phase transitions (QPT) 
in the equilibrium shape as a function of nucleon num- 
ber. This has led to the proposal [3, |4j and empirical 
verification @, Q of a new class of models, called critical 
point symmetries (CPS). These, in turn have spurred an 
abundance of experimental searches 0, (8|, Q for nuclei 
satisfying the predictions of these CPS as well as investi- 
gations into the presence of quasidynamical and partial 
dynamical symmetries at the critical point [lfj . Ill| . 

In contrast to the usual thermodynamic phase tran- 
sitions [12j . QPT occur at zero temperature as a pa- 
rameter in the Hamiltonian is varied [13j. They are 
attracting much attention in a variety of physical sys- 
tems, including Josephson-junction arrays and quantum 
Hall-effect systems [14J . The transition from Bardeen- 
Cooper-Schricffcr (BCS) pairing correlations [15| to 
Bose-Einstein condensation (BEC) as a result of in- 
creasing the strength of the pairing interaction in ul- 
tracold alkali atoms has been recently seen experimen- 
tally [lg, LL2I, LUl • In nonrigid polyatomic molecules, a 
transition from rigidly-linear to rigidly-bent shapes has 
also been studied [19j . A relationship between QPT and 
quantum entanglement has also attracted much atten- 
tion [20, [2IJ . Central to this active field is the search for 
simple empirical signatures of QPT and the identification 
of their order. The properties of QPT in atomic nuclei 
are therefore of quite broad interest, especially since nu- 
clei are finite-body systems in which the number of con- 
stituents can be varied experimentally and theoretically. 

It is the purpose of this Letter to identify a new, easy- 
to-measure, observable in one such system, atomic nuclei, 
that acts as an order parameter for QPT in nuclei and 
which can distinguish first and second order phase tran- 
sitions. Further, we will show that a particular value of 
this observable identifies the entire critical region, even 
for structures not satisfying a specific CPS. Finally we use 
this value to point to the possible existence of a hereto- 
fore undefined new symmetry at the critical point. 

In the context of an algebraic approach to nuclear 



+ CN 

O 

W 



+ -H 




FIG. 1: (Color online) The ratio E{Q+)/E{0+) as a function 
of C, for three values of \ for (a) Nb = 15 and (b) Nb = 100. 
The inset to (a) shows the corresponding behavior for v' 2 [281 ] . 



structure, the concept of phase transitions was developed 
a number of years ago using the intrinsic state formal- 
ism [1J, [23, l23j of the Interacting Boson Approximation 
(IBA) model [24j. This model, constructed in terms of 
the group U(6), has three dynamical symmetries corre- 
sponding to different nuclear shapes: a spherical nucleus 
that can vibrate [U(5)]; an ellipsoidal deformed axially 
symmetric rotor [SU(3)], and an axially asymmetric ro- 
tor [0(6)]. Nuclei exist which manifest these symme- 
tries but most nuclei deviate from them. Spherical to de- 
formed transition regions in the IBA from U(5) to SU(3) 
and U(5) to 0(6) undergo, in the large boson number 



(large valence nucleon number) limit, first and second 
order phase transitions, respectively. 

More recently, critical point transitions have been de- 
scribed in a geometrical framework, in terms of the Bohr 
Hamiltonian. These new CPS, E(5) and X(5) 0, cor- 
respond to second and first order phase transitions be- 
tween a vibrator and a rotor, differing in the 7 degree of 
freedom. In both cases, an infinite square well potential 
in the deformation, (3, allows for analytic solutions. 

The critical points in the IBA characterize the points in 
shape transitional regions where various observables [251 ] 
such as -R4/2 = E(4i )/E(2f) or electromagnetic transi- 
tion strengths, B(E2; 2f — » 0^), Q invariants [261) . as 
well as a measure of the wave function entropy [27[ vary 
most rapidly: their first derivatives have an extremum 
and their second derivative reverses sign. 

With increased interest in QPT in nuclei, several in- 
vestigations have sought order parameters to distinguish 
first and second order phase transitions. Two of these 
are [28| ^2, the difference between the expectation val- 
ues of the number of d-bosons in the IBA, n^, in the 
first excited + state and the ground state, and v' 2 , pro- 
portional to the isomer shift between the first 2 + state 
and the ground state. Both have sharp changes in phase 
transitional regions and exhibit a different behavior for 
first and second order phase transitions (see inset to Fig. 
1(a)) for small Nb- 

However, these differences do not persist in the large 
Nb thermodynamic limit [29j |. In a more practical vein, 
Ref. |30| considered the ratio of electromagnetic transi- 
tion strengths, B i/2 = B(E2;4:+ -> 2+)/B(E2;2+ -v 0+), 
in the IBA which exhibits a peak prior to the critical re- 
gion for the U(5)-SU(3) transition, similar to v' 2 . Ref. 
[311 ] determined that a difference between U(5)-SU(3) 
and U(5)-0(6) transitions persists for -B4/2 and another 
B(E2) ratio up to Nb = 50. These signatures however, 
are often difficult to measure, and have substantial un- 
certainties, particularly far off stability. Moreover, the 
B4/2 ratio varies only slightly as a function of structure, 
typically in the range 1.4-1.7 only. 

The present work presents four key results. We identify 
an order parameter based on the energies of the first ex- 
cited + state and the 6 + level of the ground state band. 
We show that the ratio E(6f)/ 'E(0 2 ) can distinguish be- 
tween first and second order phase transitions and that 
this distinction persists in the large Nb limit. We also 
identify, for the first time, a signature of the entire criti- 
cal region which closely characterizes the line of the first 
order phase transition in the IBA in the large Nb limit, 
that is a signature which identifies a nucleus as lying near 
the critical point, regardless of the 7-dependence of the 
potential. Finally, we also discuss a possible underlying 
symmetry in the critical region of the IBA. 

Calculations for this study were performed using the 
extended consistent Q [32|, (33| formalism (ECQF) of the 
IBA with a Hamiltonian given by [3J| 
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FIG. 2: (Color online)(a) Experimental E(Q+)/E{Q+) ratio 
as a function of neutron number for the Nd, Sm, Gd, and Dy 
isotopes, (b) Same for the Xe and Ba isotopes. For smaller 
neutron numbers, the Of state was taken in the ratio if its 
B(E2) decay was consistent with the a — N — 2 state. This 
corresponds to N = 74 in Xe and N = 76,78 in Ba. 
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where hd = eft ■ d, Q x 
the number of valence bosons, and c is a scaling factor. 
The above Hamiltonian contains two parameters, £ and 
X, with C, ranging from to 1, and x ranging from to 
— v7/2. In this parameterization, the three dynamical 
symmetries are given by £ = 0, any x f° r U(5), £ = 1, 
X = -a/7/2 for SU(3) and C = 1-0, X = 0.0 for 0(6). 
Calculations were performed with the code IBAR [35] 
which allows boson numbers up to 250. 

In Fig. 1, the results of calculations for the ratio 
£(6+)/£(0+) are given for N B = 15 and 100. The calcu- 
lations include the U(5)-SU(3) transition (x = — v7/2), 
the U(5)-0(6) transition (x = 0.0) and an intermediate x 
value. This ratio is 1.5 in the U(5) limit and approaches 
zero in the large Nb limit of SU(3) as the 2 state rises in 
energy. For the case of the U(5)-0(6) transition (x = 0), 
we present two sets of calculations. For Nb = 15, above 
C ~ 0.6, the dotted line follows the + state belonging to 
the a = N-2 family of the 0(6) symmetry whereas, the 
dash-dotted line follows the + belonging to the a = N 
family (the actual 2 level). 
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FIG. 3: (Color online) (Top) Line of degeneracy between the 
0+ and 6+ levels (solid red) for N B = 10, 40, 100, and 250 
in the IBA triangle. (Bottom) Line of degeneracy between 
the + and 10+ levels (solid blue) for iV fl = 250 (left) and 
between the + and 14 + levels (solid green line) for Nb — 
250 (right) in the IBA triangle. The dashed lines denote the 
critical region in the IBA obtained in the large Nb limit from 
the intrinsic state formalism 13, 22, 231. 



For small N Bl Fig. 1(a), E(6^)/E(0%) exhibits a mod- 
est peak before the first order phase transition point (% = 
— Vy/2, solid curve) followed by a sharp decrease across 
the phase transition. For the second order case, it maxi- 
mizes at U(5) and gradually decreases with increasing £. 
This behavior is identical to that for v' 2 in Fig. 1(a) [in- 
set]. The intermediate x value also exhibits a small rise 
before the phase transition. Thus, as x approaches the 
second order case, the unique features of the first order 
phase transition diminish, but are still present, although 
they may be difficult to distinguish in finite nuclei. 

Figure 1 (b) illustrates that the above behavior persists 
and is enhanced in the large Nb limit. The features of 
_E(6^)/£'(0j) are almost identical to the behavior of the 
B4/2 ratio discussed in Ref. [30| . The transition region 
becomes sharper and occurs for a narrower range of £ val- 
ues, while the peak prior to the phase transition increases 
in magnitude. With increasing Nb, the -E(6^)/-E(0^) ra- 
tio has values close to the £ = and ( = 1 limiting cases 
for ( values just outside of the transition region. 

It is also important to determine if this behavior is ob- 
served in nuclei. In Fig. 2 (a), the S(6 + )/-E(0j) ratio 
for the Nd-Dy isotopes is plotted as a function of neu- 
tron number. The Nd, Sm, and Gd isotopic chains are 
well established as undergoing a transition from spher- 
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FIG. 4: (Color online) Energies of low-lying states (nor- 
malized to 25(2 + ) = l) of the Hamiltonian of Eq. (1) with 
X=-\/7/2, (=0.4729, an d N B =250. The parameter ( was 
chosen to reproduce the approximate degeneracy of -E(0 + ) 
and E{6+). 



ical to deformed shapes, with 150 Nd, 152 Sirn and 154 Gd 
lying close to the phase transitional point [6|,LZl,l!| and ap- 
proximately manifesting X(5). Interestingly, as in Fig. 1, 
E(6i )/.E(0j) peaks just prior to the phase transitional 
nuclei, followed by a sharp drop occurring for the pro- 
posed critical point nuclei. (In 150 Nd, the rise occurs 
somewhat earlier.) The Dy isotopes behave similarly to 
the Sm and Gd isotopes, although data on N < 86 is 
lacking. Figure 2 (b) shows the E(6^)/E(0^3) ratio for 
the Xe and Ba isotopes which lie in a 7-soft region be- 
tween the U(5) and 0(6) symmetries that contains only 
a second order phase transition. Consistent with this, no 
discontinuity is observed. 

One key signature of the X(5) critical point symme- 
try is a near degeneracy of the 6 + level of the ground 
state band and the first excited + state. Close in- 
spection of Fig. 1(b) shows that the i?(6+)/S(0j) ratio 
passes through a value of 1.0 close to the phase transition 
point. This suggests that a degeneracy between E(6^) 
and -E(Oj) is related to the critical point. We note that 
this is reflected empirically in Fig. 2, where the N = 
90 nuclei also have .E(6 + )/.E(0 + ) very close to unity. In 
Fig. 3(top), the line oiE(0£) = E(6+) is given in the IBA 
triangle for Nb — 10, 40, 100, and 250. For each boson 
number, there is a single trajectory in the triangle which 
corresponds to this degeneracy condition. With increas- 
ing Nb , this line quickly approaches the phase transition 
region of the IBA, intersecting it almost exactly for Nb 
= 250. Thus, it appears that a degeneracy of -E(O^) and 
E(6f ) is a signature of the line of first order phase tran- 
sitions in the IBA in the large Nb limit. 

Degeneracies are often associated with underlying sym- 
metries. It is therefore of considerable interest that the 
critical point of the IBA in the large Nb limit exhibits 
the same degeneracy between -E(Oj) and -E(6+) as ap- 
proximately given by X(5). To further investigate the 
presence of underlying symmetries in the phase transition 
region of the IBA, Fig. 4 gives a detailed level scheme 
from the IBA with x = ~v7/2 and C = 0.4729 (which 



approximately reproduces the degeneracy of .5(0+) and 
.E(6+ )). This point lies very close to the critical point, 
Ccrit = 0.4721, obtained in the exact infinite Nb limit. 
Clearly, additional degeneracies are also present. Succes- 
sive ground band members with J > 2 with J/2 odd, 
are nearly degenerate with successive excited + states 
[e.g., (6+, 0+), (10+ 0+), (14+ 0+)]. Included in Fig. 3 
(bottom) are lines of degeneracy between the 0+ and 10+ 
levels (left) and the 0+ and 14+ levels (right). They too 
persist along the entire phase transition region. 

The continuation of the degeneracies, that appear at 
the intersection of the critical line and the U(5)-SU(3) 
leg of the triangle, into the interior of the tri ang le is rem- 
iniscent of the quasi-dynamical symmetries [3a, [33] asso- 
ciated with the U(5)-0(6) and U(5)-SU(3) transitions in 
the IBA, as well as in the vibrator to 7-unstable rotor 
transition in the Bohr model [38| . If the degeneracies in 
the interior can indeed be explained by quasi-dynamical 
symmetries, it would significantly expand the relevance 
of such ideas. In the present case, of course, the under- 
lying symmetry, exact or approximate, has not yet been 
identified and needs further investigation. 

In conclusion, we have identified a simple, empirical 
observable, the ratio £'(6+)/-E ! (0+) of the energies of the 
6+ state and the 0+ state which can serve as an order 



parameter identifying phase transitional behavior, and 
whose behavior can distinguish between first and second 
order phase transitions. Experimental data on the Sm 
and Gd isotopes, which are well known to exhibit first 
order phase transitional behavior, exhibit exactly the be- 
havior predicted by the IBA for the i?(6+)/i?(0+) ratio. 
Nuclei in second order transitional regions do not. 

Beyond its applicability as an effective order parame- 
ter, it was also found that a particular value of the ratio, 
S(6+)/-E(0+) = 1.0, provides a signature of the critical 
region which persists along the line of a first order phase 
transition cutting across the triangle in the IBA in the 
large Nb limit. The predicted degeneracy of -E(6+) and 
i?(0+) in the IBA critical region for large Nb is consis- 
tent with the X(5) critical point symmetry. Finally, the 
possibility of an underlying quasi-dynamical symmetry, 
reflected in a recurring pattern of degeneracies, in the 
critical region of the IBA was discussed. 
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